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Abstract 

In this paper some useful formulas are developed to evaluate integrals 
having a singularity of the form (t-x) _m , m>l. Interpreting the integrals 
with strong singularities in Hadamard sense, the results are used to obtain 
approximate solutions of singular integral equations. A mixed boundary value 
problem from the theory of elasticity is considered as an example. Parti- 
cularly for integral equations where the kernel contains, in addition to the 
dominant term (t-x)" m , terms which become unbounded at the end points, the 
present technique appears to be extremely effective to obtain rapidly con- 
verging numerical results. 

1. Introduction 

The mixed boundary value problems in physics and engineering may gener- 
ally be expressed in terms of a "singular" integral equation of the form 
» 

k f ( JL) = 9(2L) > 2L e0 (1) 

. 

D 

where g is a known bounded function and the kernel k is usually singular. 

The nature of singularity of k is dependent on the choice of the density 
function f in formulating the problem. For example, in one dimensional 
integral equations arising from potential theory if f(t) is selected to be 
a "flux", then k has an ordinary Cauchy singularity (t-x)~^. On the other 
hand if f is a potential, then k has a strong singularity of the form 
(t-x) . Particularly in two dimensional integral equations, formulating 

the problem in terms of a potential rather than a flux type quantity has 
certain advantages. Because of this it is worthwhile to develop effective 
techniques for evaluating singular integrals with strong singularities. In 
actual physical problems the density function f is either bounded or may 
have integrable singularities on the boundary of D. Thus, in one dimensional 
integral equations the integral on the left hand side of (1) may be inter- 
preted in Cauchy principal value sense for a Cauchy kernel, whereas it would 
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be unbounded in the case of a strong singularity (t-x) . Despite this 

in the latter case the physical problem can still be solved provided the 
integral is interpreted in Hadamard sense by retaining the finite part 
only. 


The concept of finite part integrals was first introduced by Hadamard 

[1] in connection with divergent integrals of the form 
x 


f(t) 

(x-t) p+i 


dt , 


( 2 ) 


where f is bounded and p is an integer. In spite of this relatively early 
beginning, the adoption of the concept in applications has been rather 
slow [2,3] . It is mainly due to Kutt's work [4]-[6] that in recent years 
the idea is finding relatively wide applications. To demonstrate Hadamard's 
basic idea we consider the following integral 


S 0 (x) 


b 

— ^-r = 2(b-x)* , (x<b) , 
j x 


(3) 


from which, differentiating both sides separately, it follows that 

b 


dx S o (x) " 2 


dt 


1 


(t-x) T 2 (t-xY- 


1 


t=x 


(b-xV 


— > (x<b) . 


(4) 


In (4) it is seen that the derivative of S Q (which is bounded) is the dif- 
ference between a divergent integral and an unbounded integrated term. 
Noting that the integrated term is independent of b, we may now consider 
the derivative of S Q as being the "finite part" of the divergent integral 
and define 


I 

I 


dt 

(t-x) 


t 2 


lim [ 
c->x 


dt 


(t-x)* (c-x) : 


X] = - 


( t>— x) ■ 


(x<c<b) 


(5) 


Following are some other examples : 
b b 


1 dt - +1 = lim [ 

K (t-x) a c->x 


dt 


1 1 


(t-x) 


ct+1 


a (c-x) a ' 


:3 = 


- ± (b-x) 


-a 


(a>0), (6) 
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"f FT = lim [f FT + lo g( c “ x )] = log(b-x) , 

l_x c-*x 
x c 


L 

f 


dt 


(t-x)‘ 


F^x " x^a * (a<x<b) ’ 


_d_ 

dx 


b 

f(t)log|t-x|dt = - | dt > (a<x<b) 


a 

b 


i \ dt - f * , (a<x<b, 


a 

b 


4- I — ^T+T dt = { 1 T C— ! -s+T-Mt)* , (a<x<b , a>0) 
^ T (t-x)“ +1 L (t-x) a+1 


(7) 

(8) 


(9) 


( 10 ) 


01 ) 


In this paper first some useful formulas for the evaluation of certain 
singular integrals are developed. The results are then used to obtain 
effective numerical solutions to integral equations having kernels with 
strong singularities and some examples are given. 


2. Evaluation of Finite Part Integrals 

With an eye on applications to one dimensional mixed boundary value 

problems, in this section we will describe some simple techniques for evalu- 

-2 

ating the finite part integrals having (t-x) as the kernel. Let F(t) be 
a bounded function with continuous first and second derivatives and the 
interval be normalized such that -l<(x,t)<l. The singular integral may then 
be expressed as 


1 1 


f F(t)w(t) dt . 

f [F(t)-F(x)-(t-x)F' (x)] 



Il (t-x) 2 

ii 

(t-x) 2 


1 

i 



+ F(x)| 

willdt + F , (X )1 wCtldt 
Ct-x) 2 ^ t_x 

(-1<X<1) , 

(12) 
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where w(t) is the fundamental function of the corresponding mixed boundary 
value problem and may be determined by using a suitable function theoretic 
method [7]. For simple physical problems w is given by 

w(t) = 1, (l-t 2 ) +i , (l-t) H . (13) 


One may now note that the first integral on the right hand side of (12) is 
bounded (the integrand approaches F"(x)w(x) as t+x) and the remaining 
integrals may be evaluated by using the following expressions: 


-1 

1 

I 

h 


-1 


dt _ 
t-x 



dt 

(t-x) 2 


1 1 
1-x " 1+x ’ 


VTTt 7 

t-x 


dt = -ttX, 


(-1 <X<1 ) , 


1 

f 

-1 

1 

* 

« 

-1 


y^t 7 

(t-x) 2 


dt 


-7T 


dt 

(t-x)/] 7 ! 7 


= 0 


(-1 <x<l ) , 


(-1 <X<1 ) , 


(14) 

(15) 

(16) 

(17) 

(18) 


-1 


dt 

(t-x) 2 /l-t z 


= 0 , (-1<x<1) , 


f w dt = log B) ’ (x<1) 


-1 

1 


(x<1) , 


(19) 


( 20 ) 


( 21 ) 
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1 

f dt = 1 

^1 (t-x) /1-t /1-x 

1 


log B , ( x<l ) , 


tt^= = t=x ( " ra + lv / r5‘ 109 B) » (x<1) » 


J, (t-xp/Pt 


( 22 ) 


(23) 


I /Pt 

J- (t-x) : 


_ 2/^ (-l) m , 2m-5 

m /_ t\/-, ,.\ ft 2 (m- 1 ) ( 1 


(m-1) (l-x)(l+x) 


=rf 

-1 


** dt, 
m-1 ’ 


(t-x) 


dt 




(x<l) , 
1 


(24) 


J ml (t-x) m /Pt 


✓?(-!)" + 2m- 3 I dt 

(m-1 )(l-x)(l+x) nv "^ 2 (m-1 ) ( 1 -x) ft _ x) m-l ’ 


(x<l) , 


(25) 


where m is an integer (m>2) and 

B = 0 + y ^)/0 . 


(26) 


In solving integral equations it is often convenient to express the 
unknown function F(t) in terms of a polynomial with undetermined coefficients. 
In such problems the following expressions may be quite useful: 


I 

f W*" 2< 5„< x > > 


(27) 


f dt - - ZJ r^ &Q„w-<W*n ■ 


-1 

1 


(t-x) ‘ 


1-x 


^ dt = T n+1 (x) , (n>0) , 


(28) 


(29) 
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( 30 ) 


r ] u n (t)/nt^ 

J — (t-x) z — dt = ’ 7T ^ n+1 ^ u n ( x ) ’ ( n >°) • 

- 1 T n (t) JO, (n=0) 

(t-x)/l-t z l TrU n _-j ( x) , (n>l) , 


(31) 


T n (t) 


f 


0 , (n=0,l) 


I 

f, (t-x) 2 4^ Qt 1 U„(x) +2J1 U n . 2 (x)], (n>2), (32) 


where P n , Q n and T n ,U n are the Legendre and Chebychev polynomials of 
first and second kind, respectively. Also 


,n 


-1 

1 


B n« - T f_ dt ' Jo b *f ’ ( "- 0) 


c n( x ) " it f ( t-x)^ dt “ k f Q c k x * ( n -°) 


-1 

n rvl - if 1 t n dt . f° • (r T 0) 
n 77 ^ (t-xJ/Pt 2 * \l Q d k x * (n - ]) 


(33) 


(34) 


(35) 


1 n (■ 0 , (n=0,l) , 


(36) 


l$(x) = =f dt = i L\) x n -^ l 

" ^ (t-x) X TTFl X m ^ 


/Pt 

(t-x)" 1 


dt 


n-X 


+ E A X x n " X " k , ( x<l , n>0) , 
k=0 K 


(37) 
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1 


t n dt 


dt 


s £w - + 


-1 


(t-x) A /Pt 


X 

s 


m=l 


/ n x n-x+m 
v x-m ; x 


1 

!, 




n-X i n-l-k 

+ E Bj x n A K , (x<l, n>0) , 
k=0 K 


( 38 ) 


where the coefficients b^, c k> d^, e^, A^, B^ and the expressions for the 
polynomials B n , C n , D n and E n for n=0,...,5 may be found in Appendix A. 

In (37) and (38) X is a positive integer, (. n ) is the binomial coefficient 
and the integrals in the summations can be obtained from (20) -(26) . 

Even though there are also Gaussian type integration formulas developed 
by Kutt [5] for the evaluation of the singular integral 
b 



they are not very convenient for solving integral equations by using the 
standard quadrature method which requires the use of fixed stations t- , 
since in Kutt's formulas t^ vary as x is changed (see also [8]). 


3. Solution of Integral Equations 

Let us now assume that the mixed boundary value problem is reduced 
to the following one dimensional integral equation: 
b 

£k s (t,x) + k(t,x)]f(t)dt = g(x) , (a<x<b) , (40) 

a 

where the kernel k is square integrable in [a,b] and g is a known bounded 
function. If the unknown function f is a "potential" type quantity, then 
the singular kernel k s has a strong singularity (i.e,, it contains terms 
of the order (t-x)“ n , n>l). The fundamental (or the weight) function w(t) 
of the problem may be determined from k s and f may be expressed as 
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f(t) = F(t)w(t) , a<t<b 


9 


(41) 



where F is an unknown bounded function. In solving the integral equations 
with strong singularities the application of quadrature formulas do not 
seem to be very practical. In these problems the simplest and the most 
effective technique appears to be to approximate the unknown function F 
by a truncated series as 

N 

F(t) = E a a (t) , (42) 

n=0 n n 


and to determine the coefficients a p by a weighted residual method. Here 
<j> n may be any convenient complete system of functions. Substituting from 
(41) and (42) into (40) we obtain 


N 

£ a G (x) = g(x) , (a<x<b) , 
n=0 n n 


where 


G„(x) 


4 k s (t»x)$ n (t)w(t)dt + 


b 

k (t»x)<{> n (t)w(t)dt . 
a 


(43) 


(44) 


The coefficients a p may then be determined from the following system of 
algebraic equations: 


2 a G n (xU.(x)w,(x)dx = 

n=0 n J n J J 

3 


b 

g(x)pj(x)wj(x)dx, (j=0,l 
a 


,N) , (45) 


where p. is a coordinate function in a complete system (e.g. , a set of 
orthonormal polynomials) and Wj is the corresponding weight. The functions 
<f> n (t) and 4ij(x) are usually selected in such a way that their orthogonality 
properties may be utilized. In practice one may use trigonometric func- 
tions, Legendre polynomials, Chebychev polynomials, delta functions or 
any linearly independent set of polynomials such as t n and x'L Quite 
clearly the numerical work in (45) may be reduced considerably if we 
select 

Wj(x) = 1, *j(x) = 5(x-xj) , (j=0,l,.,.,N) . 
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(46) 



By doing so we can use a sirrple collocation method to reduce (43) to the 
following algebraic system: 

N 

E a n G n( x i) = 9 ( x ,• ) » (j=0,l ,. . . ,N) . (47) 

n=0 n J J 

Although the collocation points x- can be selected arbitrarily, in general 

SJ 

they are chosen as the roots of Legendre or Chebychev polynomials. Even 
though there is no restriction on the choice of x- , a synrotric distribu- 

J 

tion with respect to the origin with more points concentrated near the 
ends seems to help. One may also note that in case of a resulting ill- 
conditioned system one could select (M+l) coordinate functions with 

J 

M>N in (45) or (46) and determine (N+l) unknowns a n from a set of (M+l) 

equations by using the method of least squares. 

Needless to say, if the integral equation (40) contains only a domi- 
-1 -2 

nant kernel (t-x) or (t-x)~ , one may always obtain the closed form solu- 
tion by expanding the functions g(x) and F(t) into appropriate series 
and by using the results given in the previous section and Appendix A. 

4. Application: A Crack in an Infinite Strip 

In fracture mechanics the problem of an infinite strip containing a 
crack perpendicular to its boundaries has been of wide interest since 
this geometry can be used as an approximation to a number of structural 
components and laboratory specimens. The related boundary value problem 
will be discussed below and the numerical treatment of the resulting 
singular integral equation will be given to demonstrate the solution tech- 
nique that was outlined in the previous section. 

As shown in Fig. 1, the crack lies perpendicular to the stress-free 
boundaries and is under prescribed surface tractions p(x). The problem 
requires solving the Navier's equations 

v 2 u. + U|^ = 0 , (1=1 ,2 ) (48) 

subject to 


- 9 - 



cr 1 1 ( 0 ,y ) = a 12 (0,y) = a-j^h.y) = a-| 2 (h,y) = 0 , (-«><y<oo) , (49) 


a-| 2 (x>0) = 0 > (0<x<h) , 
a 22 (x,0) = p(x), (a<x<b) 

u 9 (x,0) = 0 , (0<x<a, b<x<h) 

L ) 


(50) 


(51 a,b) 


where u-| , u 2 are the x,y components of the displacement vector, is the 
stress tensor referred to x,y coordinates and < is an elastic constant 
( k=3-4v for plane strain and k=(3-v)/(1+v) for plane stress, v being the 
Poisson's ratio. The stress and displacement components are related through 


°1j = u(u i,j +u j,i ) + XU k,k S ij ’ (1 ’ J " k=1 > 2) 


(52) 


where y and x are Lame's constants. The solution of (48) satisfying the 
symmetry condition (50) may be expressed as [9] 


u i (x,y) 


uj(x,y) = 


uf + u] + u? , (i=l ,2) , 


a 


(53) 

(54) 


u 2 (x,y) 


1 

2tt(1+ic) 


f v(t)[( K -i) £ + 


(55) 


r 2 = (x-t) 2 + y 2 , 


(56) 


00 

u](x,y) = -| t A ] + (q + x)A 2 ]e" ctX cosayda 


(57) 


■ 7 ! 


(A^+A 2 x)e _ctX sinayda 


(58) 
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u?(x,y) ■ - § 


[B-| + (j + h-x)B 2 ]e~^ h ~ X ^ coseyde 


u 2 (x,y) = \ f [B-,+B 2 (h-x)]e“ B ^ h " x ^ sineydg 


(59) 


(60) 


where V is the auxiliary function defined by 


V(x) = u 2 (x,+0) - u 2 (x,-0) , (a<x<b) 


(61) 


cl 2 

In this solution u.. , u^ and u^ are respectively associated with an infinite 
plane with a crack and the half planes x>0 and x<h. Using the homogeneous 
boundary conditions (49), the unknown functions A-j(a), A 2 (a), B-j(b) and 
B 2 (b) can be expressed in terms of V(x) and the mixed boundary conditions 
(51) may be shown to reduce to the following integral equation [9] : 



V(t)k(t ,x)dt 


-IT 



)p(x) , a<x<b , 


(62) 


where the kernel k(t,x) is given by 

k(t,x) = k^(t,x) + k^(h-t,h-x) + k 2 (t,x) + k 2 (h-t,h-x) 


k i (t ’ x) ■ ~ Tt+*F + * 

k 2 (t,x) = [f 1 (t,x,a)e"“ (t+x) + f 2 (t,x,o)e“ a ^ 2h+x-tl ]da , 


(63) 

(64) 

(65) 
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f-|(t,x,a) = ^- e 2an {8a 4 h 2 tx-12ci 3 h 2 (t+x)+2a 2 [9h 2 +h(t+x)+tx] 

- 3a[2h+t+x]+5+e 2a ^[-2a 2 tx+3a(t+x)“5j} , 

f 2 (t,x, ) = jj {-4a 3 [hx(h-t)]+6a 2 [h 2 +h(x-t)] 

+ o[-l 0h+t-x]+34€ _2ah [a(x-t) -3] } , 

D = l-(4a 2 h 2 +2)e“ 2ah +e' 4ah . (66a-x) 

Note that for h-*° vanishes and the integral equation for the half plane 
is recovered. 

Normalizing the interval (a,b) by defining 
t = (^)r + ( 6 ± 5 .) , x = (^)s + (6±S-) , 

V(t) = (^)f(r) 

the integral equation (62) becomes 

1 1 
* 

dr + f(r)K(r,s)dr = g(s), -1 <s <1 , (69) 

-1 

where 

K(r,s) - (^2.) 2 k(t,x) , g(s) = »(^)p(x) . 1 (70) 

The cases a>0 and a=0 represent the internal and the edge crack, respec- 
tively, and these two problems will be treated separately. In each case 
the solution will be assumed to be of the form 

f (r) = F(r)w(r) (71) 

where the fundamental solution w(r) can be determined from the dominant 



(67) 

( 68 ) 
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behavior of the singular kernels in the integral equation and is found 
to be 


w(r) = /l-r z 

, internal crack , 

(72) 

w(r) = /l-’r 

, edge crack. 

(73) 


Internal Crack : a>0 

Following the procedure described in the previous section, F(r) is 
now approximated in terms of a truncated series of Chebychev polynomials, 

N 

F(r) = E a U I (r) , (74) 

n=0 n n 

By substituting from (71), (72) and (74) into (69) and by using (30) we 
obtain 


N 

E a n [-7r(n+l)U n (s)+h n (s)] = g(s) , -!<s<l , 


where 


h n (s > 


1 

-1 


U n (r)K(r,s)/T 


dr 


(75) 


(76) 


The unknown coefficients a n are then determined from equation (75) by 
selecting a convenient set of collocation points such as 


WV = 0 » s j = cos(^L J.) , (j=0,l,. ,.,N) , (77) 

Once the solution is obtained, the stress intensity factors which are the 
main parameters of interest in fracture problems, can be calculated from 


k-j(a) = lim /2(a-x) c^U* 0 ) > (x<a) 
x*>a 


= (Jrr) lim 


V(t) 
t-*a /2(t-x) 


(t>a) 


= (-^-) 
V+l ; 


b-a 


F(-l) 


(78) 
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( 79 ) 


k-, (b) 


lim /2(x-b) a 2 2( x »°) > x >b 
x-»b 




v(t) 


2y % b-a 


= (-^-1 
V+l' 


t->b /2(b-t) 

F(l) 


, t<b 


Equations (78) and (79) are obtained from (62) by observing that the left- 
hand side in (62) gives the stress component (^(xjO) outside as well as 
inside the cut (a,b). 

Table 1 shows the stress intensity factors for an internal crack in 
a half-plane under uniform loading, p(x) = -Po as an example. 


Table 1. Normalized stress intensity factors for an internal 
crack in a half-plane. (N+l) terms are used in 
approximating the unknown function. 


( b+a l 

k-j ( a) 

k-j (b) 

N+l 

^b-a J 

p 

p o 2 

p 

K o 2 

1.01 

3.6387 

1.3298 

15 

1.05 

2.1547 

1.2536 

10 

' 1.1 

1.7585 

1.2108 

10 

1.2 

1.4637 

1.1626 

6 

1.3 

1.3316 

1.1331 

6 

1.4 

1.2544 

1.1123 

6 

1.5 

1.2035 

1.0967 

4 

2.0 

1.0913 

1.0539 

4 

3.0 

1.0345 

1.0246 

4 

4.0 

1.0182 

1.0141 

4 

5.0 

1.0112 

1.0092 

4 

10.0 

1.0026 

1.0024 

4 

20.0 

1.0006 

1.0006 

4 
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Edge Crack : a=0 

The solution of the integral equation (62) for a=0 needs more care. 
This is due to the fact that the kernel k(t,x) becomes singular as t and 
x approach 0 simultaneously (similarly in (69) K(r,s) becomes unbounded 
as r and s approach -1.) 

For a weight function A -r certain relations involving singular inte- 
grals of power series have been presented in Section 2. Therefore, if 
we express the unknown function F(r) as 

N 

F(r) = E a r n , (80) 

n=0 n 


the singular integrals may be evaluated from (37) by letting x=2. The 
integral equation (69) now becomes 


N 

1 a n M s ) = ’ -1<S<1 , 

n=0 n n 


where 


3 n (s) ’ f dr + 


-1 


r n /1-r K(r,s)dr , 


■1 


or using the notation of (37), 

1 


S n (s) = R n 2 (s) + 


^/T^r K(r,s)dr , 


-1 


(81) 


(82) 


(83) 


The integral in (83) can be evaluated numerically, however, as s-v-1 , the 
value of the integral becomes unbounded. It may be observed that for s=-l 
R n 2 te) is also unbounded resulting in a bounded value for G n (-1). To 
determine the coefficients a n the collocation points may be selected as 
in (77). 

For h-*» the kernel K(r,s) is simply 


K(r,s) 


. 1 . 12(s+l) 12(s+l ) 2 

(r+s+2) 2 ( r+s+2 ) 2 “ (r+s+2) 4 


(84) 
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from which it follows that 


G n (s) = R2(s)-R2(-s-2)+12(s+1)R3(-s-2)-12(s+1) 2 R^(-s-2) . 
In the limiting case we find 

G n (-l+0) = (4n+l )(-l ) n . 



(85) 


( 86 ) 


(87) 


As a first example we again consider a semi -infinite plane with an 
edge crack. In this case the kernel of the integral equation is given in 
closed form (see (69) and (84)) and the numerical analysis can be carried 
out quite accurately. For a uniform crack surface pressure p(x) = -P 
and for various values of N the calculated stress intensity factor k(b) 
and the relative crack opening displacement V(0) are given in Table. 2. The 
table also shows the correct value of k(b) which was calculated from the 
infinite integral given in [10] (see Appendix B). It is seen that the 
convergence of the method is extremely good. 

* 

The second example is concerned with a long strip of finite width h 
which contains an edge crack of length b and is subjected to a uniform 
tension P Q (p(x)=-P 0 ) (table 3) or pure bending M (p(x) = - (1 - ^-) 

(table 4) away from the crack region. In the numerical analysis the num- 
ber of collocation points was increased until the accuracy of the last 
significant digits given in tables 3 and 4 were verified. In no case 
more than 20 points were needed. 

Aside from providing accurate answers to some very practical ques- 
tions, the results given in Tables 3 and 4 are important in that they follow 
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Table 2. Normalized stress intensity factor and crack open- 
ing displacement for an edge crack in a half-plane 


N+l 

k-j(b) 

V(0) 

P 0 /F 

M+k' PqD 

1 

1.062652 

1.502816 

2 

1.126950 

1.423476 

3 

1.124283 

1.457747 

4 

1.121818 

1.455918 

5 

1.121442 

1.454520 

6 

1.121451 

1.454224 

7 

1.121483 

1.454211 

8 

1.121504 

1.454241 

9 

1.121514 

1.454264 

10 

1.121518 

1.454278 

15 

1.121522 

1.454298 

20 

1.121522 

1.454298 


1.121522* 



*The correct value of stress intensity factor (calculated 
from the infinite integral given by Koiter [10] . 
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Table 3. Normalized stress intensity factor and crack open- 
ing displacement for an edge crack in a strip under 
uniform tension. 


b/h 

¥ b) 


o 

a. 

-> 0 

1.12152226 


0.00001 

1.121522 

0. 14543xl0~ 4 

0.001 

1.121531 

0.1 4543x1 0 - ^ 

0.1 

1.1892 

0.15490 

0.2 

1.3673 

0.36543 

0.3 

1.6599 

0.70358 

0.4 

2.1114 

1.3048 

0.5 

2.8246 

2.4702 

0.6 

4.0332 

4.9746 

0.7 

6.3569 

11.246 

0.8 

11.955 | 

31.840 

0.85 

18.628 

63.288 

0.9 

34.633 

158.94 

0.95 

99.14 

708.8 
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Table 4. Normalized stress intensity factor and crack open- 
ing displacement for an edge crack in a strip 
under pure bending, a =6M/h 2 


b/h 

k-j(.b) 

( 2u i V(0) 

a vHT 
n 

'1+K h 

+ 0 

1.12152226 

0.14543xl0" 4 

0.00001 

1.1215 

0.001 

1.1202 

0.14535x10”^ 

0.1 

1.0472 

0.14529 

0.2 

1.0553 

0.31822 

0.3 

1.1241 

0.56141 

0.4 

1.2606 

0.94130 

0.5 

1.4972 

1.5924 

0.6 

1.9140 

2.8387 

0.7 

2.7252 

5.6432 

0.8 

4.6764 

13.989 

0.85 

6.9817 

25.990 

0.9 

12.462 

60.965 

0.95 

34.31 

253.7 
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routinely from the technique presented in this paper, for very deep cracks 
(t»0.8h) are not available in literature, and are extremely difficult to 
obtain by using other methods. For example, the solution of the corres- 
ponding singular integral equation having a Cauchy type dominant kernel 
by using a Gaussian integration formula requires much greater computational 
effort than the technique presented here for the same accuracy and for 
b>0.8h has an extremely slow convergence. 
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APPENDIX A 


The coefficients given in the equations (33) - (38) : 


0 , for n-k = even. 


b, = 


l , r(^) 


2/7 r(Hiil) 


(n>0, n>k) , 
, for n-k = odd , 


c k 


0 , for n-k = odd 

(k+1 ) 2 } 

1 ST ZW6 » ' for n-k 

2vV r(^P) 


(n>0, n>k) , 


= even 


d k* 


0 , for n-k = even 


= < 


r 

1 , r(n^) 


(n>l , n>k) 


/7r(^4il) 


, for n-k = odd 


<■ 


^ 0 , for n-k = odd 

l (k+n f , k . 

*V r(^) ’ ° r " eve " 


(n>2, n>k) 


^ c-D J (j) ^ • OeP) 


B k = 2^ <-D J (j> 2j+T ’ ("1°) 

where (") is the binomial coefficient and for the gamma function we have 

r(-±) = -2r(i) = -2/F. 

Some examples of the polynomials B n , C n , D R and E n are given below 
[see (33)-(36)]: 
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-X 


9 


C o = -> 


9 


= -x 2 + 


= -2x , 


-x 3 + ^ 
a 2 9 

_ X 4 + xi + 1 

X 2 8 » 

_ x 5 + xi + X 

X 2 8 ’ 

_ X 6 + x^L + xi + J_ 

X 2 8 16 » 


- -3x 2 + j , 
= -4x 3 + x , 


= x 2 + ± , 

= X 3 + f , 

= X 4 + xi + 3 
X 2 8 * 


= 3X 2 + -7T 9 


= 4x 3 + x . 
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APPENDIX B 


Stress Intensity Factor for the Edge Crack 
Calculated from Koiter's Results 


The edge crack problem in a semi -infinite plane has been considered in 
the literature many times and mostly for comparison reasons. The stress 
intensity factor 1.1215 has become a standard when comparing numerical 
techniques for the solution of singular integral equations. For uniform 
pressure p 0 applied on the crack surface, a closed form expression for 
the stress intensity factor in terms of an infinite integral is given by 
Koiter [10]: 



p 0 vfr Jn A 


where A is calculated from 


1 og( “ ■ si nh 3 2 ) da 

v'B z +a z [coshira-2a 2 -l ] 


and B is an arbitrary constant greater than 1. 

The result is independent of the choice for B and numerical calcula- 
tions show that 


— — = 1.12152226 , 

Po^ 

where there may be an error only in the last digit. 
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